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The brief description of installation TAU-2 intended for long-term monitoring of the half-life value 
T (T 1 / 2 ) of the ^^^Po is presented. The methods of measurement and processing of collected data are 
reported. The results of analysis of time series values r with different time step are presented. Total 
of measurement time was equal to 590 days. Averaged value of the half-life was obtained 

r = 163.46 ± 0.04 yis. The annual variation with an amplitude A = (8.9 ± 2.3) • 10”"^, solar-daily 
variation with an amplitude Aso = (7.5 ± 1.2) • 10“^, lunar-daily variation with an amplitude 
Al = (6.9 ± 2.0) • 10”“^ and sidereal-daily variation with an amplitude As = (7.2 ± 1.2) • 10”“^ 
were found in a series of r values. The maximal values of amplitude are observed at the moments 
when the projections of the installation Earth location velocity vectors toward the source of possible 
variation achieve its maximal magnitudes. 

PACS numbers: 27.80.-l-w, 23.60.-|-e 


I. INTRODUCTION 

At the last time in works intended to search for limits 
of the realization of the decay constant conservation law, 
a level of sensitivity not less than 2 • 10“"^ was reached for 
several radioactive isotopes. In the work [l[ the authors 
showed an amplitude of a possible annual variation of 
the ^®®Au half-life {T 1/2 = 2.69445 days), that was mea¬ 
sured with the relative uncertainty of ±7 • 10“® does not 
exceeds ±2 • 10“"* of the central value. Variations with 
periods from several hours up to one year were excluded 
at the level of 9.6 ■ 10“® (95% C.L.) during the measure¬ 
ments of the ^^^Cs half-life {T 1/2 = 10942 days) in the 
0. The annual variation was excluded at the level of 
8.5 • 10“^ (95% C.L.). Variations of an activity with pe¬ 
riods of 3-150 days were excluded at the level of 2.6 ■ 10“^ 
(99.7% G.L.) during the measurement of the activity 
in the Q . It was shown that an amplitude of the annual 
variation does not exceeds of 6.1 ■ 10“® (95% C.L.). Vari¬ 
ations of an activity with periods less then one year were 
excluded at the level of 4 • 10“^ during the measurement 
of the ^^^Th activity in the 

A count rate of the detector recording the source ra¬ 
diation was a subject of investigations in the all men¬ 
tioned works. A high sensitivity of the measurements 
was reached by using of a relatively high count rate 
(~ 10^ s~^); of a control and a stabilization of conditions 
of the measurements and by using of additional arrange¬ 
ments for shield of the set-ups from outer background. 

At the same time, the evidences of a presence of the 
annual variations of different effects caused by a radia¬ 
tion of the investigated isotope are cited in a series of 
articles. For example, characteristics of the annual vari¬ 
ations of the count rates of the detectors used to many 
years measurements of the ^^Si and ^^®Ra sources activ¬ 
ities discussed in the work Q. The amplitudes of varia¬ 
tions are equal to ^ 1 • 10“^. The authors have examined 


possibilities of an appearing of such variations as a result 
of seasonal variations of the detectors characteristics or 
the one of an annual modulation of the isotopes decay 
rates themselves under the action of an unknown factor 
depending of the Earth-Sun distance. Results of continu¬ 
ous measurements of the decay rates of the ^°®Ag, ^^^Ba, 
i^^Cs, i^^Eu, ®5Kr, 226Ra and ®°Sr sources made 

in the Physikalisch-Technische Bundesanstalt (PTB) was 
discussed in the works d, Q • Statistically significant an¬ 
nular variations with the amplitudes of (6.8 — 8.8) • 10“"* 
were observed in the data for all this isotopes. The au¬ 
thors made a comparison between spectral power func¬ 
tions obtained from the data and the observed radial os¬ 
cillations of the Sans surface. However, the researchers 
from the PTB point in the [3] at the laboratory factors ca¬ 
pable to cause similar variations and warn against hasty 
conclusions about possible new physical phenomena. 

It is clear that any conclusions about a possible new 
physical effect could be made after the complete exclu¬ 
sion of variations caused by influence of the known ter¬ 
restrial geophysical, climatic and meteorological factors 
on the source-detector couple count rate. Not all such 
factors could be detected and took into account during 
the measurement and data processing. For example, an 
annual variation with the amplitude of (4.5 ± 0.8) • 10“^ 
was found as a result of a processing of the data col¬ 
lected at 500 days of Earths surface measurement with 
source in the Q. It was found that this variation 
corresponds completely to the known annual variation of 
the cosmic rays intensity and could be explained by a 
cosmic rays background events contribution to the total 
detectors count rate. A variation with the ^ 300 days 
period and 4 • 10“® amplitude was found in the data col¬ 
lected at 480 days in the underground measurement with 
the source. It was found that this variation corre¬ 

late with a variation of a daily averaged dead time per 
event and could be explained by a modulation of the RG 
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circuit providing the shaping time of the amplifier. 

The weak point of the experiments aimed to monitor 
a stability of a controlled radiation count rate is the high 
sensitivity to the similar variations of the measurement 
conditions. It seems that this shortcoming became unim¬ 
portant for the decay constant determination based on 
direct registration a life time of nucleus between a birth 
and a decay. This methodic was realized by us in the 
Q for the ^^^Po which decays with 164.3 fis half-life @ 
by emitting the 7.687 MeV a-particle. This isotope ap¬ 
pears mainly in the exited state 87%) in the 
/3-decay. Half-lives of the exited levels does not exceed 
0.2 ps [l 3 | and they discharge instantly with regard to the 
scale of the half-life. Energies of the most intensive 
7 -lines are equal to 609.3 keV (46.1% per decay), 1120 
keV (15.0%) and 1765 keV (15.9%). So, the /3-particle 
and 7 -quantum are emitted at the moment of a birth of 
^^"‘Po nuclear (start) and the a-particle are emitted at 
the decay moment (stop). Measurement of ’’start-stop” 
time intervals allows one to construct decay curve at an 
observation time and to determine the half-life from its 
shape. The ^^®Ra source {T 1/2 = 1600 years) was used 
as a generator of nuclei which arise in the decay 

sequence of the mother isotope. 

The direct measurement of a nuclear life time allows 
one moreover to study the radioactive decay law itself. 
The theoretical models discussed in the [ll|, [Tl| predict 
that the decay curves could deviate from the exponential 
law in the short- and very long-time regions of the time 
scale. The theoretically predicted |13l - [lq so called quan¬ 
tum Zeno effect consisting in a slowing down of the decay 
rate in a case of constant observations at the decaying 
object presents a special interest. Experimentally Zeno 
effect was proved [l^ in repeatedly measured two-level 
system undergoing Rabi transitions, but not observed in 
spontaneous decays. 

A limitation of the annual variation amplitude was 
set at the level of 3.3 ■ 10“^ at the first stage of our 
measurements Q. Factors limiting a sensitivity were 
revealed and ways of its optimization were designed. 
The improvements of the set-ups, methods of measure¬ 
ments and data processing, results of an analysis of pos¬ 
sible sources of systematic errors were described in the 
work m. The annular variation with an amplitude of 
(6.9 ± 3.0) ■ 10“^ and the solar-daily variation with an 
amplitude of ( 10.0 ± 2 . 6 ) • lO”'^ were obtained as a result 
of a processing of the data collected at 480 days of the 
measurements. 

A brief description of the used installation and results 
of a broadened analysis of the data collected at 590 days 
are given in the present work. 

II. METHOD OF MEASUREMENTS 

The TAU-2 set-up consists of the two scintillation de¬ 
tectors D1 and D2 used in the work. The D1 was made 
of two glued discs of a plastic scintillator (PS) with the 


18 mm diameter (d) and 0.8 mm thickness (h). A thin 
transparent circular bag glued of two 2.5 fim lavsan lay¬ 
ers placed between the discs. A radium spot is deposited 
preliminary in the center of inner surface of a one lavsan 
circle. The D1 registers ,5-particles from the ^^“^Bi de¬ 
cays and a-particles from the ^^^Po decays. The massive 
detector D2 made of Nal(Tl) crystals intended for the 
7 -quanta detection. Two Nal(Tl) crystals (d=150 mm, 
h=150 mm) placed by ends one to another with the 10 
mm gap. The D1 is placed into a gap between D2a and 
D2b. The light collection is fulfilled from a lateral side of 
the PS disc installed into deep narrow well with a reflect¬ 
ing wall. The measurements are carried out in the low 
background room of the underground laboratory DULB- 
4900 of the BNO INR RAS at the depth of 4900 m of 
water equivalent in the additional shield made from 
Pb (15 cm). 

A registration of the pulses in the set-up is carried 
out by the two-channel digital oscilloscope La-n20-12PCI 
which is inserted into a personal computer (PC). Pulses 
are digitized with 6.25 MHz frequency (160 ns/channel). 
The DO pulse recording starts by a signal from the D2 
which registered ^^“^Bi decays 7 -quanta. A D2 signal 
opens a record of a sequence with 655.36 /iS total du¬ 
ration where first 81.92 /rs time is a ’’prehistory” and the 
last 573.44 /xs is a ’’history”. Duration of a ’’history” 
exceeds the three ^^^Po half-lives. 

A scintillation detector D1 in the TAU-2 has a rela¬ 
tive a//3 light output ^ 0.1 [l^. As result the pulses 
from the a- and /3-particles have the comparable am¬ 
plitudes. This circumstance was used to a preliminary 
selection of the ’’useful” events by the ’’on-line” program 
prepared the data for a PC recording. It is obtained a 
number of pulses in the D1 and D2 channels and its de¬ 
lays. The pulses have tested at correspondence to the 
’’correct event” criterions such as 1 ) the only one pulse 
in the D2 channel, 2) presence of the prompt coincid¬ 
ing pulse in the D1 channel, 3) presence of the only one 
delayed pulse in the ’’history” of the D1 channel and 4) 
absence of any pulses in the ’’prehistory” of the Dl. The 
frames do not correspond to the correct event criterions 
are rejected. A ’’correct” event appearing time, pulse 
amplitudes and its appearing time are recorded into the 
PC memory. This information allows one to ’’off line” 
process the data for the different regions of pulse am¬ 
plitudes. A count rate of the ’’right” events is equal to 
^ 12 s“^. A rate of an information accumulation is equal 
to 25 Mb-day“^. 


III. SEARCH FOR LONG DURATION 
VARIATIONS OF THE TAU-2 DATA 

The spectra of the /3-pulses (spectrum 1) and a-pulses 
(spectrum 2) of the Dl detector and of the 7 -pulses (spec¬ 
trum 3) of the D2 detector corresponding to the selection 
conditions mentioned above are shown in FiglT] The peak 
at the channel ~ 1600 on the spectrum 2 is formed by the 
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FIG. 1. Spectra of the coincided in 573.44 jis delay of pulses 
from the detector D1 (/3-spectrum 1) and D2 ( 7 -spectrum 3) 
and the a-spectrum 2 of one of the first pulse D1 to install 
TAU-2. 



FIG. 2. The start-stop delays distribution obtained in 590 
days at the TAU-2. The decay-curve of the ^^^Po with r = 
163.45 ± 0.04 fis. 



FIG. 3. Dependence in time of the r-value obtained at the 
TAU-2 with the week step. 



Time, week 

FIG. 4. Distribution in time of the ^^^Po half-life values ob¬ 
tained by the ’’moving-average method”. Gurve 1 approxi¬ 
mation by function @ with the parameters: A = 5.4 • 10~^, 
uj = 27r/365 d“^, = 83 d (since the 1st of January). 

Gurve 2 approximation by function © with the parameters: 
A = 8.9 • 10“'‘, uj = 27r/365 d~\ ^ = 174 d. 


7.69 MeV a-particles. A total time of the data collection 
is equal to 730 days at the period October 2012 October 
2014. A decay curve constructed for the total data set 
is shown in Fig© A value of the r obtained from this 
data is equal to 163.45 ± 0.04 /rs. The r-value equal to 
(164.58 ±0.29 (stat.)±0.10 (syst.)) us was measured at 
the Gran Sasso in a recent work [^. These two values 
are compatible in the Icr limits. 

The continuous in time data set was divided at por¬ 
tions of the equal length to search for possible time vari¬ 
ations. The decay curve has constructed for the each 
partition data set and the continuous in time sequence of 
the r-values with the specified time step has found. The 
time dependence of the r-values with the week time step 
is shown in Fig© The r-values were defined for the 3.2- 
560 IIS delay-time region by means of a y^-approximation 
of the decay-curves collected at one week each with an 


exponential function 

y = a ■ exp{—ln{2) ■ t/r) + b. (1) 

The moving-average method (moving summation) was 
used to search for a possible time variations of the r- 
values. A time interval with the duration equal to about 
0.5 of the expected period is choose to search for any 
harmonic component and the r-value is determine for this 
interval. Than, the interval have shifted for a one step 
and the procedure is repeated. The 0.5 year interval and 
the one week step were chosen to search for the annular 
variations. A decay curve was constructed for the data 
collected at 0.5 year and the r-value was determined for 
it. Than, the interval have shifted for the 1 week and a 
determination was repeated. The result of the analysis 
is presented on the Fig© 
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FIG. 5. Distribution in time of the ^^^Po half-life values for 
one week data sets. Curve 1 approximation by function m 
with the parameters: A = 8.9-10 = 27r/365 d ^,(p = 174 

d. Curve 2 - time dependence of E.-S. distance (right scale). 

Harmonic component have presented in the data as it 
seen from the figure. The approximation dependence 

T{t) = To ■ [1 + • sin{oj ■ {t + (^))] (2) 

where A = 5A- 10“"^, w = 27r/365 d“^, tp = 83 d (since 
the 1st of January) obtained from the data by means of 
X^-method is shown by red color on the Fig0] too. It 
is easy to show that the original dependence of the week 
data set has the same period (1 year), the amplitude large 
at 7r/2 time and shifted at 0.25 year (0.5 of the sliding 
interval) than the approximation one. The annual wave 

r(t) = To • [1 -I- 8.9 • 10“^ • sin{u! ■ {t + 174))] 

obtained by such a way is shown on Fig (S] together with 
the week data set of the r-values. 

It is follows from the consideration that the periodic 
component with the 1 year period and the amplitude of 

A = (8.9 ±2.3) ■ 10“'‘ 

is presented in the data. The maximum of the periodic 
function is observed at the ^October 22. Phases of the 
obtained period and the period of the Earth-Sun distance 
variation differ at 3 months. Therefore, obtained periodic 
variation could not be explained by the distance chang¬ 
ing. However, a parameter of the Earth velocity relative 
to the Sun is exists in connection with the Earth-Sun 
(E.-S.) distance annular variation. A dependence of this 
parameter is shown of Figl5] (curve 2). The obtained an¬ 
nular wave of ^^^Po half-life (FiglSl curve 1) is coincides 
in phases with the E.-S. velocity dependence within the 
±1 week accuracy. 

The described above method was used to search for 
the daily variation in the solar, lunar and sidereal times. 
A days duration was divided on 24 hours in the each 
case. The durations of the lunar and sidereal days in 
the solar time are equal to 24 h 50 min 28.2 s and 23 h 



half-life values for a week data sets. Curve 1 approximation 
by function © with the parameters: A = 8.9 -10 cu — 
27r/365 d“^, ip = 174 d. Curve 2 - time dependence of E. to 
S. velocity. 



FIG. 7. Solar-daily variation of half-life of ^^^Po obtained by 
means of moving-average method. The solid curve - approx¬ 
imation of function ([2]) with the parameters: A = 4.8 • 10“^, 
ui = 27r/24 h“^, p = —6 h. 

56 min 4.09 s correspondingly [2l|. The 12 h duration 
period was chosen as an averaging interval. All events 
registered in the 0-12 hour interval at the total period 
of the measurements were used to construct the decay 
curve. A half-life value determined from this curve. The 
interval was moved for the 1 hour and procedure was 
repeated after that. The result of a search for the daily 
variation in the solar time is shown on the Figl7]with the 
parameters {A = 4.8 • 10“^, w = 27r/24 h“^, ip = —6 h) 
of the approximation of function © . The daily variation 
is described good enough by this sine function as it seen 
from the Fig IT) A reconstructed sought for dependence 

r(t) = To • [1 ± 7.5 • lO”"^ • sw((27r/24) • t)] 

is shown on the Fig[5] (curve 1) together with a depen¬ 
dence of the Earth surface point velocity relative to the 
Sun due to the Earth rotation (curve 2). Their phases 
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FIG. 8 . Curve 1 - Solar-daily variation of half-life of 
approximated by function m with the parameters: A — 7.5 • 
10“"^, u) = 27r/24 h“^, (^ = 0 d. Curve 2 - Earth surface point 
velocity relative to the Sun due to the Earth rotation. 


are coincide in the range of ±0.5 hour. The amplitude of 
a solar daily variation is equal to Aso = (7.5 ± 1.2) • 10“'^. 

The result of a search for the daily variation t(<) in 
the lunar time is shown on the Figl^ (curve 1) with the 
parameters of the approximation function Q: 

A = 4.4 • 10-^ w = 2^24 h-\ if = 6 hu, 

where h-id lunar day hour. The reconstructed sought for 
dependence 

r(t) = rO • [1 ± 6.9 • 10"^ • sm((27r/24) • (t ± 12))] 

is shown on the FiglU (curve 2) together with a depen¬ 
dence of the Earth surface point velocity relative to the 
Moon due to the Earth rotation (curve 3). Phases of the 
curves (2) and (3) are coincide in the range of ±0.5 hour. 
The amplitude of lunar daily variation is equal to 

Al = (6.9 ±2.0) • 10"^ 

It was done a test of a real presence of the lunar 
daily variation in the data set. A starting point of the 
data summation was shifted at 15 days. A phase of 
the variation should shifts at ^ 12 hours because of a 
difference in a time duration of a solar and lunar days 
[(50 min 28.2 s)T5 Ri 12 hj. The result is shown on 
FigHni A phase had shifted at 12 hours as it was waited 
for a real lunar daily variation. 

The result of a search for the daily variation r(t) in 
the sidereal time is shown on the Fig llll (curve 1) with 
the parameters of the approximation function Q : 

A = 4.6 • 10-^ w = 2^24 h"/, ^=-6 hia, 

where hsd sidereal day hour. The reconstructed sought 
for dependence 


r(t) = rO • [1 ± 7.2 • lO”"* • sm((27r/24) ■ t)] (3) 


FIG. 9. Lunar-daily variation of half-life of ^^“^Po obtained by 
means of moving-average method. The data set starting point 
is 24:00, 30 September 2012. Dash dot curve - approximation 
function eq. m with the parameters {A — 4.4 ■ 10 ■*, ca = 
2txI2A 1/3 = 6 hid). Dash curve - sought Lunar-daily 
variation r{t) = vo • [1 -|- 6.9 ■ 10“"^ • sm((27r/24) ■ {t -I- 12))]. 
Solid curve - Earth surface point velocity (Ve.p.m/) relative 
to the Moon due to the Earth rotation. 
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FIG. 10. Lunar-daily variation of half-life of ^^“^Po obtained 
by means of moving-average method. Starting point shifted 
for 15 days (15 October 2012). 


is shown on the FigfTTl (curve 2) too. The amplitude of 
a sidereal daily variation is equal to 

As = (7.2 ± 1.2) • 10“^. 


It was done as test of a real presence of the sidereal 
daily variation in the data set. A starting point of the 
data summation was shifted at 182 days. A phase of 
the variation should shifts at ~ 12 hours because of dif¬ 
ference in a time duration of solar and sidereal days 
[(3 min 55.9 s)-182 rs 12 h]. The result is shown on 
Fig HD A phase had shifted at 12 hours as it was waited 
for a real sidereal daily variation. 
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FIG. 11. Sidereal-daily variation of half-life of ^^“^Po ob¬ 
tained by means of moving-average method. Curve 1 - func¬ 
tion eq. m with the parameters {A = 4.6 -10 a; = 27r/24 
hJd 1 ‘P = —& had). Curve 2 - sought SD variation r(t) = 
TO • [1 + 7.2 • lO""* • sm((27r/24) • t)]. 


FIG. 13. Anti-Sidereal-daily variation of half-life of ^^^Po 
obtained by means of moving-average method, (dasd = d + At 
= 24 h 3 min 55.9 s). 


IV. RESULTS AND DISCUSSION 
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FIG. 12. Sidereal-daily variation of half-life of ^^“^Po obtained 
by means of moving-average method. Starting point shifted 
for 182 days. 


The possibility of a stochastic realization of sidereal 
variation was tested also. An analysis was repeated for 
an artificial anti-sidereal day. Its duration in the solar 
time is equal to 24 h 3 min 55.9 s. The result of the anal¬ 
ysis is shown on the Fig ITS] Any variation does not seen 
within the statistical errors. It gives a confidence that 
the sidereal daily variation is present really in the data. 
It is difficult to correlate the sidereal daily variation with 
any cosmic object at present time because of a nature of 
the possible influence is unknown. Nevertheless, if such 
object is exists than it should be an annular variation 
connected with it. The question needs an additional fur¬ 
ther investigation. 


The moving-average method was used to search for an¬ 
nular and daily variations because of a need to improve 
a statistic of analyzed decay curves and to increase mul¬ 
tiply a sensitivity of the analysis to a value of a possible 
variation. As it seen from a comparison the data on Fig|3] 
and 01 a value of a statistical error was decreased at 5 
times by increasing of a data accumulation interval from 
the one week up to 26 weeks (0.5 year). Using of the 
Fourier-method and wavelet-method to the analysis is 
complicated by a relatively short size of the time series. 
Besides, a time r-series could include not only annular 
variation but semiannual, for example. The integrated 
dependence lost the unambiguity in this case. The used 
method targeted at a search for the specific variation 
is determined a univocal binding exactly. It is neces¬ 
sary to compare obtained set of the r-values with similar 
data measured independently on the other similar set-up 
to find a confidence that the variations are nonrandom. 
Such measurements are carried out in the BNO INR RAS 
with the TAU-1 set-up since May of 2014. The set-up is 
placed in the low background laboratory ’’CAPRIZ” at 
the 1000 m w.e. depth. The statistics needed to carry¬ 
ing out of a comparison of coincidence behavior of the 
T- values time sets is accumulates. It is clear that an es¬ 
sential improvement of statistics could be achieved by a 
multiple increasing of the dataset rate. The quadratic in¬ 
crease of a random coincidence fraction will occur for the 
used ^^'*Po source in a case of a considerable growth of its 
activity. The value of random coincidence is equal to 1% 
at 12 s“^ dataset rate. Such large value connected with a 
high total activity of all ^^®Ra daughter elements in the 
source and relatively long half-life of ^^"'^Po. Because of 
it, an increase of dataset rate for the ^^"'^Po without rela¬ 
tive growth of a random coincidence background is possi¬ 
ble by means of increasing of the number of independent 
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measuring set-ups. This variant seems as a hardly feasi¬ 
ble. Another possibility could be realized by using of pair 
of radioactive isotopes with a similar decay scheme but a 
shorter half-life. The pair ^^^Bi~^^^Po (ri /2 = 4.2 /is) is 
a good candidate for such source. This isotopes are the 
daughter products of ^^®Th (ri /2 = 7340 y) [T^ which 
would be used as a generator isotope. The preparation 
of the new set-up for ^^®Th is carry out at the BNO INR 
RAS at present time. 


V. CONCLUSIONS 

The results of analysis of the data obtained with TAU- 
2 set-up at the new step of measurements are shown on 
the presented work. The set-up is intended to carrying 
out of a long duration control at a value the ^^'^Po half-life 
constant. The methods of measurement and processing 


of collected data are reported. Results of the analysis of 
time series values of t with different time step are pre¬ 
sented. Total time of measurements was equal to 730 
days. Averaged at 590 days value the ^^"^Po half-life was 
found ri/2 = 163.46 ± 0.04 p,s. It is shown that the con¬ 
stant feels the daily and annular variations of the un¬ 
known nature. The annual variation with an amplitude 
A = ( 8.9 ± 2 . 3 ) ■ 10 “"'^, solar-daily variation with an am¬ 
plitude Aso = ( 7.5 ± 1 . 2 ) • 10 “^, lunar-daily variation 
with an amplitude A^ = ( 6.9 ± 2 . 0 ) • 10 ““^ and sidereal- 
daily variation with an amplitude As = ( 7.2 ± 1 . 2 ) ■ 10 “'’’ 
were found in a series of t values. The maximal values 
of amplitude are observed at the moments when the pro¬ 
jections of the installation Earth location velocity vectors 
toward the source of possible variation achieve its maxi¬ 
mal magnitudes. The measurements are continuing. 

The work was made in accordance with INR RAS and 
V.N.Karazin KhNU plans of the Research and Develop¬ 
ments. 
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